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ABSTRACT 


This  Note  discusses  the  design  of  software  to  be  used 
for  the  planning  of  minor  planet  observations  and  the  identi¬ 
fication  of  a  particular  asteroid  when  it  is  observed  at  a 
given  time.  In  each  case  the  large  number  of  cataloged 

<,•>  -y  *  Jr  t<t 

asteroids  (<£  2200)  has  required  that  the  software  be  optimal 
with  respect  to  execution  speed.  This  is  accomplished  by 
using  an  approximate  geocentric  ecliptic  position  to  eliminate, 
as  quickly  (in  terms  of  CPU  time)  as  possible,  the  minor 
planet  from  further  consideration.  In  particular,  using  the 
type  of  analysis  applicable  to  near-stationary  artificial 
satellites,  (and  carefully  quantifying  the  nature  of  the 
approximations)  an  approximate  geocentric  ecliptic  latitude 
and  longitude  are  obtained  in  such  a  fashion  that  the  computa¬ 
tions  needed  for  these  are  not  only  minimal  but,  should  the 
actual  location  of  the  minor  planet  be  needed,  maximally 

(  £  r-  *  f 

used  there  too.  The  accuracy  is  (for  e  <_\  0.3,  i  <  30°) 

,  i 

3?9A'^)27lrin  geocentric  ecliptic  longitude  and  1?8, 

in  geocentric  ecliptic  latitude.  With  these  values  we  can 
eliminate  approximately  half  of  all  asteroids  because  they 
will  be  on  that  half  of  the  celestial  sphere  which  corresponds 
to  daytime  .(\  Half  of  the  remainder  are  then  eliminated 
because  i_hey  will  be  below  the  observer's  (usable,  i.e., 
above  altitude  30°)  horizon.  Here  the  asteroid's  zenith 
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distance  is  computed  in  the  ecliptic  coordinate  system.  The 
user's  interaction  with  the  software,  complete  documentation 
of  the  calculations,  a  user's  guide,  sample  outputs,  and  the 
rationale  of  observing  asteroids  are  also  discussed. 
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I.  INTRODUCTION 

With  the  advent^-  of  regular  asteroid  observing  at  the 

Experimental  Test  System  (hereinafter  ETS)  of  the  Ground- 

based  Electro-Optical  Deep  Space  Surveillance  (GEODSS) 

* 

program  it  became  clear  that  a  planning  aid  for  routine 

observations  and  an  identification  aid  for  both  routine 

observing  and  minor  planet  searches  were  necessary.  The 

planning  aid  (whose  acronym  became  ASTPT  from  ASTeroid 

PoinTing)  would  be  a  non-real  time  program  similar  to  the 

one  used  at  the  ETS  for  scheduling  the  observations  of 

2 

artificial  satellites  .  The  identification  aid  (ASTID) 
would  be  a  program  executed  in  real  time.  Since  most  asteroids 
are  in  essentially  the  same  orbit  (i.e.,  circular,  in  the 
plane  of  the  ecliptic,  2.7  A.U.  from  the  Sun)  the  identif ication 
space  should  be  as  large  as  possible.  Hence  provision  to 
use  position,  angular  velocity,  apparent  magnitude,  color 
index,  and  lightcurve  period  would  be  necessary.  Again 
there  is  a  similar3  real  time  program  at  the  ETS  but  it  only 
uses  position^. 

* - 

GEODSS  will  supplant  the  Baker-Nunn  camera  system. 

^Angular  velocity  information  and  in  and  out  of  orbital 
plane  information  (relative  to  the  supposed  artificial 
satellite)  is  provided  to  the  operator  but  not  automatically 
used  for  discrimination  purposes.  The  other  data  is  unavailable 
for  most  high  altitude  artificial  satellites  at  present. 


The  difficulty  is  the  large  size  of  the  asteroid  file  - 

4 

the  current  one  contains  2188  entries  and  is  rapidly  growing  . 
The  current  GEODSS  (essentially  high  altitude)  artificial 
satellite  file  has  502  entries.  Moreover,  execution  of  the 
ETS  correlation  software  frequently  requires  20-25  seconds. 
Clearly  something  better  in  concept  is  needed. 

This  Note  discusses  the  design  of  fast,  efficient  soft¬ 
ware  which  fulfills  the  requirements  of  the  planning  function. 
Execution,  exclusive  of  the  major  output  sections,  only 
requires  ^  63  seconds.  Since  only  minor  modifications  are 
necessary  to  turn  ASTPT  into  ASTID,  and  ASTID  should  execute 
at  least  ^10  times  faster  than  ASTPT  does  (cf .  SVIX) ,  it  would 
appear  that  the  requirements  of  the  real  time  identification 
function  have  also  been  met.  A  User's  Guide  to  ASTPT  is 
provided  in  the  Appendix.  Sample  outputs  are  given  in  Figs. 

1  and  2.  Before  turning  to  the  details  of  the  software,  it 
is  appropriate  to  review  the  reasons  asteroids  are  observed 
at  all  at  the  ETS. 
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II.  WHY  WE  OBSERVE  ASTEROIDS 

We  observe  asteroids  because  pt  their  many  similarities 

with  (high  altitude)  artificial  satellites  and  by  regarding 

asteroids  as  "unknowns"  we  can  thoroughly  test  many  of  our 

operating  techniques,  methods  of  data  analysis  and  data 

reduction,  initial  orbit  generation,  and  the  differential 

correction  of  orbital  elements.  For  instance  asteroids  are 

natural  satellites  of  the  Sun  shining  by  reflected  sunlight. 

Presumably  all  of  them  (there's  detailed  data  on  'v  350  light- 

curves)  spin  about  a  fixed  axis  producing  variations  in  their 

lightcurves  which  require  analysis  for  the  determination  of 

(i)  the  sidereal  rotation  period,  (ii)  any  precessional  or 

nutational  period  that  might  be  present,  (iii)  the  orientation 

of  the  spin  axis,  (iv)  the  elimination  of  aspect  variations 

and  (v)  size  and  shape  determinations,  and  (vi)  the  elucidation 

.  5 

of  surface  mottling  .  Multicolor  photometry  (from  reference  6, 
which  is  highly  recommended  as  a  general  introduction  and 
reference  volume,  we  find  that  the  mean  B-V  for  the  744  asteroids 
with  data  is  01!1??  with  a  standard  deviation  about  the  mean  of 
O^OO)  has  proven  useful  for  separating  asteroids  into  classes 
based  on  chemical  composition.  Polarization  data  also  aids 
in  some  studies  and  the  combination  of  visible  and  infrared 
photometry  (coupled  with  assumptions  concerning  albedo  and 
emissivity)  yields  diameter  estimates.  Few  asteroids  are 


large  enough  or  close  enough  to  have  a  resolvable  disc. 

Finally,  in  this  category,  analysis  of  minor  planet  phase 

functions  may  also  yield  surface  compositions  and  an  under- 

* 

standing  of  the  opposition  effect  -  something  akin  to  this 
having  already  been  observed  for  artificial  satellites7. 

The  classical  methods  of  obtaining  positions  of  minor 
planets  have  already  been  successfully  adopted  for  real  time 
applications  at  the  ETS  {reference  8  contains  an  overview 
and  complete  references  to  supporting  documentation;  the 
current  accuracy  in  the  total  positional  error  is  =:  2 " 5 ) . 

Minor  planet  orbits  also  exhibit  strong  perturbations  due  to 
third-body  forces. 

When  we  turn  to  searches  we  find  two,  modern,  intensive, 

o 

short-term  efforts  (the  Yerkes-McDonald  Survey  and  the 
Palomar-Leiden  Survey10) .  Longer  term,  less  intensive  efforts 
are  principally  performed  at  the  Crimean  Astrophysical  Observa¬ 
tory11  (with  a  40  cm  double  astrograph)  and  by  Helin  and 
12 

Shoemaker  for  Earth-crossing  asteroids  (using  the  18" 

Schmidt  on  Mount  Palomar) .  The  results  and  techniques  of 
these  efforts  when  coupled  with  our  own  searches  for  artificial 
satellites11  and  asteroids1  should  continue  to  improve  the 
artificial  satellite  searches  in  several  areas.  Examples 

a 

The  opposition  effect  refers  to  the  abnormal  brightening 
(by  ^  0.4)  of  an  asteroid  within  a  few  (^  8)  degrees  of 
opposition . 


(and  some  successes)  include  a  one  telescope  versus  a 

coordinated  two  telescope  search13,  dead-reckoning  methods*3,14 

and  identification  procedures13.  It  is  in  the  field  of  initial 

orbit  construction  (e.g.,  the  classical  method  of  three 

15 

positions  versus  statistical  methods  ,  methods  based  on 

preliminary  distance  estimation1^,  or  methods  especially 

.  17 

adapted  to  a  particular  orbital  type  ) that  we  expect  the 
greatest  dividends. 

Clearly  the  interaction  between  the  two  sets  of  objects 
has  already  been  fruitful  and  there  is  every  indication  that 
continuing  to  mine  the  store  of  asteroid  observing,  data 
reduction,  and  compiled  information  will  further  the  goals  of 
the  GEODSS  program.  There  are  two  other  important  points. 

Since  no  other  observatory  is  as  well  equipped  as  the  ETS 
is  for  performing  the  real  time  detection  and  discrimination 
of  minor  planets,  or  for  obtaining  their  positions  in  real 
time,  the  telescope  time  is  well  spent  efficiently.  Clearly 
whatever  data  we  do  acquire  will  benefit  the  astronomical 
community.  The  other  point  is  that  no  other  observatory  is 
capable  of  performing  a  real  time  search  for  asteroids. 

Here,  especially  for  the  increasingly  important  Earth-crossing 
objects,  the  ETS  can  play  a  unique  role. 
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III.  NOTATION 


Since  asteroids  revolve  about  the  Sun  their  orbital 
elements  are  clearly  heliocentric  -  as  are  the  locations 

derived  from  them.  We,  however,  require  geocentric 

* 

(theoretically  topocentric,  see  §  VE)  positions  .  Further¬ 
more  the  ecliptic  of  1950.0  is  the  fundamental  reference 
plane  for  the  orbital  element  sets.  Therefore,  a  variety  of 
coordinates  (rectangular,  spherical;  ecliptic,  equatorial, 
and  horizon;  heliocentric,  geocentric,  and  topocentric;  1950.0 
epoch  and  epoch  of  date)  are  needed.  In  an  attempt  to  have 
the  notation  imply  as  much  information  as  possible  the 
following  scheme  has  been  adopted:  rectangular  ecliptic 
coordinates  will  be  denoted  by  x,  y,  z  and  spherical 
ecliptic  coordinates  by  r,  X,  $  (distance,  ecliptic  or 
celestial  longitude,  and  ecliptic  or  celestial  latitude). 

Their  interrelationship  is  given  by 


r  = 


r£  (X ,  $ ) 


Rectangular  equatorial  coordinates  will  be  denoted  by  x,  y ,  z 
and  spherical  equatorial  coordinates  by  r,  a,  6  (the  same 
distance  for  the  same  origin,  right  ascension,  and  declination) . 


Heliocentric  coordinates  are  symbolized  with  lower  case 
letters  (as  in  r,  a,  6  or  r,  X,  6)  and  the  corresponding 
geocentric  quantities  by  upper  case  letters  (e.g.,  R,  A,  A 
or  R,  A,  B) .  Topocentric  quantities  and  horizon  system 
coordinates  are,  it  turns  out,  unnecessary  in  practice  except 
for  the  zenith  distance  Z  (necessarily  of  date) .  Quantities 
relative  to  the  1950.0  ecliptic  or  relative  to  the  mean 
equator  and  equinox  of  1950.0  are  denoted  by  a  zero  subscript 
(e.g.,  iQ,  6q,  Aq  etc.).  The  same  quantity  relative  to  the 
mean  ecliptic  of  date  or  relative  to  the  mean  equator  and 
equinox  of  date  has  no  subscript  (e.g.,  i,  o,  A).  Distances 
are  unaffected  by  general  precession  (r  =  rQ,  R  =  Rq,  etc.). 

The  parameters  of  the  asteroid  also  have  no  subscript 
while  those  of  the  Sun  carry  €>.  Other  important  quantities 
include  the  asteroid's  absolute  B  magnitude  B(l,  0),  the 
phase  angle  f  (the  angle  in  the  Sun-asteroid-Earth  triangle 
with  vertex  at  the  asteroid) ,  the  local  mean  sidereal  time  t , 
the  observer's  astronomical  latitude  <p,  the  coefficient  of  f 
in  the  phase  function  in  magnitudes  F,  and  the  ecliptic  longi¬ 
tude  and  latitude  of  the  astronomical  zenith  (geocentric, 

epoch  of  date  only)  L  ,  B  . 

z  z 

This  system  of  notation  is  rigidly  adhered  to  for  both 
our  sakes.  Thus,  AQ  can  only  mean  the  geocentric,  ecliptic 
longitude  of  the  Sun  relative  to  the  mean  equator  and  equinox 


7 


of  date.  The  principal  symbols  are  defined  below.  When  the 
same  symbol  has  been  used  twice  the  most  important  meaning  is 
given  first. 


a  =  semi-major  axis  of  the  asteroid  orbit 
a  =  general  precession  in  longitude 
A  =  see  Eq.  (7b) ;  epoch  of  date 
A^  =  see  Eq.  (7b);  epoch  of  1950.0 
A  =  see  Eq.  (7a);  epoch  of  date 
A^  =  see  Eq.  (7a);  epoch  of  1950.0 

B(1,0)  =  absolute  blue  (B)  magnitude  of  the  asteroid 
B(rR,f)  =  apparent  exo-atmospheric  blue  (B)  magnitude 
of  the  asteroid;  see  Eq.  (15) 

Bz  -  geocentric  latitude  of  the  astronomical  zenith, 
of  date 

e  =  eccentricity  of  the  asteroid  orbit 
e^  =  lower  eccentricity  limit 
eu  =  upper  eccentricity  limit 
E  -  eccentric  anomaly;  see  Eq.  (1) 

f  =  phase  angle;  see  Eq.  (17)  and  Pig.  3 
F  =  o’I'oZS/deg;  see  Eq.  (15) 

i  =  inclination  of  the  asteroid  orbit  relative  to  the 
ecliptic,  of  date 
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i  =  inclination  of  the  asteroid  orbit  relative  to  the 
o 

ecliptic,  1950.0 
i^  =  lower  inclination  limit 
iu  =  upper  inclination  limit 

T 

Mu,v)  =  (cosvcosu,cosvsinu,  smv) 

L  =  geocentric  ecliptic  longitude  of  the  astronomical 
zenith,  of  date 

m  =  speed  of  general  precession  in  right  ascension 

m^  =  lower  apparent  magnitude  limit 

m^  =  upper  apparent  magnitude  limit 

m20  =  aPParent  S-20  magnitude;  see  Eq.  (19) 

M  =  mean  anomaly  at  time  t 

n  =  mean  motion  =  2ir/a3^2 

n  =  speed  of  general  precession  in  declination 

P  =  general  precession  matrix;  see  Eq.  (13b) 

q  =  aphelion  distance  of  the  asteroid  =  a(l  -  e) 

r  =  r  =4=*=  heliocentric  distance  of  the  asteroid 
o  o 

T 

r  =  ri(A,g)  =  (x,y,z)  *=  heliocentric,  ecliptic  location 

of  the  asteroid,  of  date 
T 

4  =  n!L{ a, 6)  =  (x,j/,z)  =  heliocentric,  equatorial  location 

of  the  asteroid,  of  date 
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R0  =  geocentric  distance  of  the  Sun 

*®  =  R0i  (AQ,B0)  =  geocentric,  ecliptic  location  of  the  Sun, 

of  date 

B©  =  R©A(A0'A©)  =  geocentric  equatorial  location  of  the  Sun, 

of  date 


(u)  =  x  -  axis  rotation  matrix  by  angle  u  = 


/I 

0 

°  \ 

0 

cosu 

sinu  I 

\o  - 

sinu 

cosu/ 

/cosv 

0 

-sinv\ 

• 

1 

° 

\sinv 

0 

cosv  / 
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R3 (w)  =  z  -  axis  rotation  matrix  by  angle  w  = 


cosw  sinw  0 
-sinw  cosw  0 
I  o  0  1 


S  =  R3  (-fiQ)R1  {-iQ)R3  (-0Jo) ;  see  Eqs.  (5) 
t  =  the  time  of  interest 

T  =  the  time  of  perihelion  passage  of  the  asteroid 


u  =  v  +  u>  =  true  argument  of  latitude,  of  date 
uQ  =  v  +  u>Q  =  true  argument  of  latitude,  1950.0 

v  =  true  anomaly,  see  Eq.  (3) 

z  =  equatorial  precession  quantity,  see  Eq.  (13c) 
Z  *  zenith  distance  of  asteroid 


a  ~  heliocentric  right  ascension  of  the  asteroid,  of  date 
aQ  *  heliocentric  right  ascension  of  the  asteroid,  1950.0 
A  =  geocentric  right  ascension  of  the  asteroid,  of  date 
Aq  «  geocentric  right  ascension  of  the  asteroid,  1950.0 
A©  *  geocentric  right  ascension  of  the  Sun,  of  date 

6  =  heliocentric  ecliptic  latitude  of  the  asteroid,  of  date 
3q  =  heliocentric  ecliptic  latitude  of  the  asteroid,  1950.0 
B  =  geocentric  ecliptic  latitude  of  the  asteroid,  of  date 
Bq  =  geocentric  ecliptic  latitude  of  the  asteroid,  1950.0 
Bg  =  geocentric  ecliptic  latitude  of  the  Sun,  of  date 
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see  Eq.  (21c) 

heliocentric  declination  of  the  asteroid,  of  date 
heliocentric  declination  of  the  asteroid,  1950.0 
geocentric  declination  of  the  asteroid,  of  date 
geocentric  declination  of  the  asteroid,  1950.0 
geocentric  declination  of  the  Sun,  of  date 
nutation  in  obliquity 
nutation  in  longitude 

obliquity  of  ecliptic,  of  date 
obliquity  of  ecliptic,  1950.0 
equatorial  precession  quantity,  see  Eq.  (13c) 
equatorial  precession  quantity,  see  Eq.  (13c) 
heliocentric  ecliptic  longitude  of  the  asteroid,  of  date 
heliocentric  ecliptic  longitude  of  the  asteroid,  1950.0 
geocentric  ecliptic  longitude  of  the  asteroid,  of  date 
geocentric  ecliptic  longitude  of  the  asteroid,  1950.0 
geocentric  ecliptic  longitude  of  the  Sun,  of  date 
see  Eq.  (2) 

local  mean  sidereal  time 

astronomical  latitude  of  the  observer 

argument  of  perihelion  of  the  asteroid  relative  to 

the  ecliptic  of  date 

argument  of  perihelion  of  the  asteroid  relative  to  the 
ecliptic  1950.0 
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longitude  of  ascending  node  of  the  asteroid  relative 
to  the  ecliptic  of  date 

longitude  of  ascending  node  of  the  asteroid  relative 
to  the  ecliptic  of  1950.0 


IV.  THE  USER'S  INTERACTION 

A.  Inputs 

While  designing  the  input/output  sections  of  ASTPT,  which 
will  primarily  be  used  to  prepare  observing  schedules,  the 
objectives  of  observing  asteroids  (§11)  at  all  have  been  the 
principal  consideration.  Thus  the  operator  can  select  apparent 
magnitude,  inclination,  and  eccentricity  ranges  and  the  order  of 
presentation  of  the  data.  For  example  one  might  only  want  to 
use  high  eccentricity  (e  >  0.23)  or  high  inclination  (i  >  14°) 
asteroids  to  test  dead-reckoning  or  initial  orbit  determination 
schemes.  To  mimic  a  search  for  near-stationary  artificial 
satellites  the  desirable  minor  planets  would  have  low  eccen¬ 
tricity  (e  <  0.23)  and  low  inclination  (i  <  14°),  and  the  out¬ 
put  should  be  arranged  in  order  of  right  ascension.  Testing 
the  efficacy  of  actual  operations  would  benefit  from  an 
angular  speed  sort.  Multi-color  photometric  data  acquisition 
and  analysis  and  light  curve  analysis  would  benefit  from  both 
apparent  magnitude  selection  and  sorting  of  the  output  by 
brightness.  Similarly  the  delineation  of  magnitude  effects 
near  the  limits  of  detection  and  discrimination  would  benefit 
from  these  options.  Therefore,  the  user  is  asked  to  specify 
six  parameters  relating  to  apparent  magnitude,  inclination,  and 
eccentricity  or  accept  their  default  values.  The  parameters 


m^  =  the  lower  (bright)  limit  for  the  apparent  magnitude 
m^  =  the  upper  (faint)  limit  for  the  apparent  magnitude 
i£  =  the  lower  limit  for  the  inclination 

iu  =  the  upper  limit  for  the  inclination 

e£  =  the  lower  limit  for  the  eccentricity 

eu  =  the  upper  limit  for  the  eccentricity 

The  default  values  are  m^  =  0m,  m^  =  18?5;  i£  =  0°,  i^ 

=  180 °;  e£  ■  0,  e  *  1.  The  bright  limit  is  arbitrary  and 
insures  that  no  minor  planet  will  be  arbitrarily  excluded. 

The  faint  limit  represents  the  experimentally  determined1 
ETS  limit  for  positive  detection  and  discrimination.  Through 
a  question  and  answer  phase  of  ASTPT  the  user  can  alter  any 
subset  of  these  six  quantities.  Illegal  entries  (e.g.,  mfa  >. 
mf,  i£  >_  iu,  e£  >  eu;  or  mf  >  18?5,  i£  or  iu  >  180°,  e£ 
or  eu  >  1;  e£  <  0,  i£  <  0)  are  prevented. 

Of  course,  the  user  must  also  supply  the  program  with 
the  date  of  interest.  The  year  is  assumed  to  be  the  current 
one  and  the  month  and  day  of  the  month  are  the  necessary 
inputs.  Specification  of  the  exact  time  of  interest  is  not 
necessary  because  the  program  automatically  produces  positions 
for  (local  mean  solar  times)  6  P.M.,  12  A.M.,  and  6  A.M. 

That  is,  if  the  user  requests  output  for  (say)  Feb.  5  he 
will  receive  data  for  6  P.M.  Feb.  5,  12  A.M.  Feb.  6,  and  6  A.M. 
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* 

Feb.  6.  If  he  insists  on  interpolating  he  can  do  so. 

Illegal  entries  are  prevented  here  too. 

After  fixing  the  date  and  making  his  choice  of  the 
apparent  magnitude,  inclination,  and  eccentricity  limits 
the  user  must  supply  the  number  of  the  tape  drive  unit  that 
the  master  asteroid  tape  file  has  been  mounted  on.  For  this 
and  all  other  operator  requested  input  the  program  immediately 
provides  an  echo  print  and  allows  the  operator  to  immediately 
modify  his  input.  Moreover  the  user,  after  entering  the  month 
and  day  of  the  month,  has  the  option  of  dealing  with  mb,  mf; 
i^ ,  iu;  e^ ,  eu;  tape  drive  unit  one  group  at  a  time  or  all  of 
them  simultaneously.  If  he  chooses  all  at  once  and  fails  to 
respond  properly  he  can  still  return  to  the  one  group  at  a 
time  mode. 

Each  cataloged  asteroid  is  assigned  a  unique  identification 
number.  Our  file  (kindly  supplied  in  machine  readable  form  by 
Prof.  Brian  Marsden,  Director  of  the  Minor  Planet  Center) 
contains  2188  entries  ordered  according  to  increasing  identi¬ 
fication  number.  This  is  also  the  standard  format  for  all 
compilations  of  asteroid  data.  Output  is  most  easily  provided 
in  this  sequence  and  represents  the  default  option  of  ASTPT. 

Even  if  other  sequences  are  desired  this  listing  will  be 

5E 

A  fast  moving  main  belt  minor  planet  would  have  a  geocentric 
angular  speed  of  l'/hr. 
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provided.  As  discussed  above  other  sequences  of  interest  are 
increasing  right  ascension,  geocentric  angular  speed,  or 
apparent  magnitude.  After  specifying  the  tape  drive  unit  the 
user  can  request  no  additional  sorting,  any  single  sort,  or 
any  combination  of  multiple  sorts  (performed  and  listed 
separately;  the  possibility  of  an  angular  speed  sort  within 
a  right  ascension  sort  may  be  added) .  Sorted  data  is  indicated 
by  the  heading  (see  Figs.  1  and  2),  not  by  rearranging  the 
output  format  to  forcefully  reflect  the  fact  that  a  sort 
had  been  performed. 

Sample  sessions  of  interaction  with  the  program  are  in 
the  Appendix.  This  includes  error  messages  too.  There  is  no 
way  of  automatically  choosing  by  default  the  default  values 
for  the  above  parameters. 

B.  Outputs 

The  form  and  content  of  the  output,  given  the  flexibility 

already  provided  to  the  user,  has  been  designed  to  minimize 

computation  (in  his  head  and  the  CPU)  and  maximize  usefulness. 

Our  earlier  experience  with  a  more  primitive  version1  of  this 

type  of  software  has  led  to  the  format  shown  in  Fig.  1.  For 
* 

each  visible  minor  planet  three  lines  are  given.  The  top  line 
has  the  geocentric  right  ascension  and  declination  (both  in 

3F 

An  asteroid  is  "visible"  if  and  only  if  it  passes  all  apparent 
magnitude,  inclination,  and  eccentricity  criteria  and  is  above 
altitude  =  30°  at  midnight. 
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sexagesimal  measure,  time  and  circular  respectively,  to  the 

*  ... 
nearest  second  )  for  the  mean  equator  and  equinox  of  midnight 

(hereinafter  the  mean  equator  and  equinox  of  date,  see  §§VB,VIA, 

VIB)  at  the  instant  of  time  corresponding  to  6  P.M.  The 

effects  of  planetary  perturbations  have  been  ignored,  the 

effects  of  parallax  [(heliocentric-geocentric  only;  not 

diurnal  (see  §VE)]  and  planetary  aberration  (§  VD)  are  included. 

The  second  line  contains  identification  number,  apparent  S- 

20  magnitude^  (including  the  effects  of  heliocentric  and 

geocentric  distance,  phase,  atmospheric  extinction,  color 

corrections  from  B,  and  the  opposition  effect;  see  §VD) ,  A 

and  A  (as  on  the  top  line)  for  the  instant  of  time  corresponding 

to  12  A.M.  ,  the  time  needed  for  the  asteroid  to  move  7 1'5  in 

minutes  of  time  to  the  nearest  minute  (which  experience  at 

the  ETS  has  shown  to  be  the  minimum  distance  for  efficient, 

positive  detection  and  discrimination) ,  the  position  angle 

(in  degrees  to  the  nearest  10°)  of  the  geocentric  angular 

velocity,  the  compass  direction  (i.e.,  NE  or  S)  of  the 

angular  velocity  to  the  nearest  22 ?5,  and  a  flag  (an  asterisk) 

* 

I  realize  that  this  is  inconsistent.  It's  also  superfluous 
(see  SIVA,  VIA)  but  it  is  the  standard  input  format  for  moving 
the  telescope  to  a  particular  position  at  the  ETS. 

f 

The  photo-sensitive  surface  of  the  ETS  cameras  is  an  S-20 
type. 
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indicating  eastward  motion  (if  appropriate)  .  The  last 
line  gives  A  and  A  (as  on  the  top  line)  for  6  A.M. 


Flags  for  high  eccentricity  or  high  inclination  were 
deemed  unnecessary.  Anyway  this  is  rare.  From  reference 
18,  for  a  sample  of  2829^  asteroids,  we  find  that  the  average 
eccentricity  in  the  asteroid  belt  is  0.153  with  a  standard 
deviation  about  the  mean  of  0.079  and  that  the  average  inclina¬ 
tion  in  the  asteroid  belt  is  7?93  with  a  standard  deviation 
about  the  mean  of  5?74.  Hence,  the  one-sigma  upper  limits 
(the  0.23  and  14°  mentioned  above)  are  accidently  comparable  to 
those  in  the  definition"^  of  a  near-stationary  artificial 
satellite.  Of  these  2829  minor  planets  only  10  had  inclina¬ 
tions  in  excess  of  30°  and  only  45  had  eccentricities  in  excess 
of  0.3.  Moreover,  high  inclination  and  high  eccentricity  seem 

i 

to  go  together  . 

Figure  2  shows  a  portion  of  a  right  ascension  sort. 

All  sorts  are  based  on  midnight  values  of  the  relevant 
parameter (s) . 

a - 

Eastward  moving  asteroids  are  in  their  retrograde  loop  and 
present  some  initial  orbit  and  dead  reckoning  problems  not 
faced  when  observing  artificial  satellites. 

+This  sample  consisted  of  a  union  of  the  (then  -1969)  numbered 
asteroids  and  the  "first-class"  orbits  from  the  Palomar-Leiden 
Survey 

•^Seventeen  minor  planets  with  e  >  0.3  also  have  i  >  22°. 
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Fig.  2.  One  page  of  ASTP  output  sorted  according  to 
right  ascension  at  midnight. 


V.  THE  COMPUTATIONS 

Some  of  the  computations  involved  are  not  executed  in 
the  rigorous  fashion.  In  this  Section  the  rigorous  (or  very 
nearly  rigorous)  procedures  are  described.  As  this  is  fairly 
complicated  (because  of  the  general  precession-parallax  inter¬ 
play  and  the  number  of  quantities  we  must  calculate)  the 
approximations  are  documented  in  §VI. 

A.  Heliocentric  Locations 

The  osculating  orbital  element  set  for  epoch  1950.0  in 

the  classical  form  (e.g.,  a,  e,  iQ,  u)q,  and  T)  is  the 

h 

starting  point.  The  common  date  of  osculation  is  0n  Dec.  27, 

1980  E.T.  As  these  are  not  perfect  and  planetary  perturbations 

are  ignored  I  expect  a  positional  accuracy  of  1'  within  a  year 

or  two  of  the  date  of  osculation.  This  is  small  enough  to 

*  t 

satisfy  both  the  planning  and  identification  requirements. 

To  obtain  the  rectangular,  ecliptic,  heliocentric  coordinates  for 
the  equator  and  equinox  of  1950.0  (r^  we  first  determine  the 
mean  anomaly  M  from 

M  =  n(t  -  T),  n  =  2ir/a 3//2 


* - 

The  field  of  view  of  the  31"  telescope  at  the  ETS  when  electron¬ 
ically  zoomed  in  the  Cassegrain  configuration  is  ^  20'. 

+If  we  view  the  asteroid  belt  as  occupying  30°  of  ecliptic  lati¬ 
tude,  then  the  mean  density  of  cataloged  minor  planets  is  0.2/sq 
degree.  Hence,  their  mean  separation  is  1?25  and  the  probability 
of  two  being  within  2'  of  each  other,  given  that  one  is  already 
there,  is  0.07%. 
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where  a  is  in  A.U.,  n  in  rad/yr,  and  with  t  (the  time  of 
interest,  i.e.,  6  P.M.  Feb.  5,  1980  local  mean  solar  time)  and 
T  being  measured  in  years  of  ephemeris  time  measure.  This 
requires  that  the  difference  between  Universal  Time  and 
Ephemeris  Time  be  known.  This  is  the  case,  in  advance,  to 
^  0?2.  Again  using  l"/min  as  representative  of  a  fast  main 
belt  asteroid  (=>  a  =  1.8  A.U.  for  a  circular  orbit)  the 
resultant  error  is  negligible. 

Next  we  solve  Kepler ' s  equation  for  the  eccentric  anomaly 
E, 

E  -  esinE  =  M  (1) 

and  obtain  the  vector 


The  computation  of  the  true  anomaly  v  and  the  heliocentric 
distance  r  is  unnecessary; 

tan (v/2)  =  [(1  +  e)/(l  -  e) ] 1/2 tan (E/2)  (3) 

r  =  a(l  -  ecosE)  (4) 

Finally  we  perform  the  rotations  to  correct  for  non-zero 
values  of  wo,  and  iQ, 
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^  =  Sp 


where 


S  =  R3(-fi0)R1(-i0)R3 <““o}  (5b) 

To  obtain  rectangular,  equatorial,  heliocentric  coordinates 
for  the  equator  and  equinox  of  1950.0  (4^)  we  need  an  additional 
rotation  by  the  obliquity  of  the  ecliptic, 


=  Rl("eo)Ho; 


=  23  °26 ' 44  '.'836 


There  is  a  very  old  trick  in  astronomy  to  simplify  the  applica¬ 
tion  of  Eqs.  (5,6).  Define  vectors  A  ,  8  and  A.  by  (with 
Su  =  sinu,  Cu  =  cosu  Vu  to  save  space) 


Coj  Cft  -  Su  Sft  Ci^  \ 

o  o  o  o  o  \ 

=  a  (Cu  Sft  +  Su  Cft  Ci  )Ce  -  Su„Si  Se„ 

oo  ooo  o  oooJ 

(Cu  Sft  +  Su  Cft  Ci  )Se  +  Su  Si  Ce  / 

\  °  °  ooo  o  ooo  J 

I  -Su  Cft  -  Cu  Sft  Ci 

o  o  ooo 

=  a(l  -  e2) 1/2  (-Su  Sft  +  CuCftCi) Ce  -  CunSi  Sen 

o  o  o  o  o  o  ooo 

(-Su  Sft  +  Cu  Cft  Ci  )Se^  +  Cu  Si  Ce^ 

o  n  n  n  n  n  n 


=  A-'c  -o' 

O 


B  =8 
-o  -o 1 


eo  *  0 


Then, 


Eo  =  <VEo'V  -  <CE  -  e>*o  +  (SE)Eo 
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io  =  Uo'  V  zo)T  *  <CE  -  e)*o  +  (SE,Eo 

As  it's  *  we  want,  and  and  8^  can  be  stored  on  the  master 
asteroid  tape,  Eq.  (8)  represents  an  efficient  method  of 
obtaining 

B.  General  Precession  and  Parallax 

If  you  introduce  an  intermediate  epoch  ("in  between"  the 

fundamental  epoch  of  1950.0  and  the  epoch  of  date;  you'll  see 

why,  and  why  in  between  is  in  quotes,  in  a  minute)  to  perform 

the  general  precession  calculation  and  start  to  enumerate  all 

the  possibilities  (rectangular  versus  spherical;  ecliptic 

versus  equatorial)  you'll  very  rapidly  wind  up  with  a  long  list 

because  it's  the  precessed  geocentric  location  we  want.  It 

turns  out  that  equatorial  is  better  than  ecliptic  (the  reasons 

aren't  given  here)  and  that  rectangular  is  much  better  than 

spherical  (no  expensive  trigonometric  functions  and  their 

equally  expensive  inverses  to  perform).  Even  so,  since  we 

could've  precessed  w  ,  ft  ,  and  i  there  are  still  five  different 

o  o  o 

methods  of  handling  the  interplay  between  general  precession 
and  the  heliocentric-geocentric  parallax. 

it 

1.  We  could  precess  wQ,fi0,  and  iQ  to  the  epoch  of  date  , 

compute  A  and  then  add  R^  to  form  R  (see  Fig.  3).  Since  we 
* 

These  formulas  are  complicated  and  can  be  found  on  pages 
273-275  of  reference  19. 
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don't  know  the  epoch  of  date  all  of  the  trigonometric  computa¬ 
tions  implicit  in  Eqs.  (5  and  6)  have  to  be  performed  within 
ASTPT .  Therefore  this  possibility  is  rejected.  For  reference, 
the  approximate  (second-order)  formulas  to  do  this  (my  advice 
is  to  go  to  possibility  #2)  are 

ft  *  Qq  +  a  -  bsin  (Qq  +  c)coti 

oj  *  u>  +  bsin(Q  +  c)csci  (9) 

o  o 

i  =  i  +  bcos  (ft  +  c) 
o  o 

where  a  is  the  general  precession  in  longitude,  b  is  the  inclina¬ 
tion  of  the  ecliptic  of  date  to  that  of  1950.0,  and  c  =  180° 

-  II  +  a/2  where  II  is  the  ascending  node  of  the  ecliptic  of  date 
on  the  1950.0  ecliptic  referred  to  the  mean  equator  and  equinox 
of  date.  This  notation  is  standard.  Numerical  values  for  a 
and  b  are 


a  =  p(t  -  1950.0),  b  =  it  (t  -  1950.0) 

where  p  is  the  speed  of  the  general  precession  in  longitude 
and  it  is  the  rate  of  rotation  of  the  ecliptic  (both  per  year 
here).  Their  approximate  values,  and  values  of  comparable 
accuracy  for  II  and  e,  are 

p  =  (50  i'256  +  01'0222T)/yr 


7T  =  (0T471  -  0?0007T)/yr 


n  =  173°57 ! 06  +  54 ! 77T 

e  =  23°27,08'.,26  -  46'.’845T 

where  T  is  the  time  since  1900.0  in  tropical  centuries.  To 

keep  Eqs.  (9)  second-order  T  should  be  (t  +  1950. 0)/2. 

2.  Precess  a)  ,  ft  ,  and  i  to  1980.0  (for  calendar  year 
o  o  o 

1980  pointing)  and  store  these  as  A'  and  8*  on  the  master 
asteroid  tape.  Now  only  (my  notation  scheme  is  running  into 
trouble) 


K*  =  (cosE  -  e)A*  +  (sinE)B'  (10) 

is  executed  in  ASTPT.  Add  to  this  R'  to  form  R'.  We  then 

— ©  — 

obtain  R  by  precessing  from  1980.0  to  the  mean  equator  and 
equinox  of  date.  Since  this  time  interval  is  short  we  can 
do  this,  with  adequate  accuracy,  via 


K A,  A)  =  £(A',  A’)  +  df  (A* ,  A')/dt  •  At 


(ID 


where,  if  t(K' ,  A')  =  (€,n,5/,  then 


(12) 


and  m  and  n  are  the  speeds  of  general  precession  in  right 
ascension  and  declination. 
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m  a  (467085  +  070279T)/yr 
n  =  ( 2 0 1* 0 4 7  -  0'.'0085T)/yr 

T  is  again  the  time  in  tropical  centuries  since  1900.0  I(t  + 

1980. 0)/2  here]  and  At  in  Eq.  (11)  is  t  -  1980.0.  This  method 
is  as  efficient  as  possibility  #3  and  these  are  the  two  most 
efficient  of  the  five  possibilities  considered  here.  The 
difficulty  is  that  R*  is  not  easily  available. 

3.  What  is  available  is  R‘  from  J an.  1,  1980  to 
Jul.  1,  1980  relative  to  the  mean  equator  and  equinox  of 
1980.0  and  R'  relative  to  the  mean  equator  and  equinox  of 
1981.0  for  Jul.  1,  1980  to  Dec.  31,  1980  (remember  for  1980 
read  any  year).  Hence,  if  we  separate  the  year  into  two 
parts ,  store  oj  ,  SI ,  and  i  in  the  A '  ,  8 '  format  on  our  tape 
for  both  1980.0  and  1981.0,  then  we  can  use  the  technique  of 
possibility  #2.  But,  when  I  say  that  R'  is  available  it  is  so 
only  in  tabular  (but  machine  readable)  form.  For  now  I  preferred 
a  simple  analytic  expression  for  R’  or  R  even  at  the  loss  of 
some  precision. 

4.  If  we  invert  the  logic  of  possibilities  2  and  3  we 
proceed  via  Eq.  (8)  for  then  add  R^  (only  available  in 
[machine  readable]  tabular  form)  to  obtain  R^  and  precess  via 

£  =  (13a) 

The  matrix  P  is  given  by 
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(13b) 


P  =  R3(-z)R2(0)R3(-Co)  =  R3  (-z-90°)R1  (0)R3  (90°-Co) 

where  £q1  z,  and  0  are  the  equatorial  precessional  quantities; 

C  =  2304  '.'95177At  +  0"30217(At)2  +  0 r01800 (At ) 3 

z  =  5  +  Oi'79303  (At)2  +  0  L'00032  (At )  3  (13c) 

0  =  2004  ”25826At  -  0 *.'42689  (At)  2  -  0 '.'04180  (At )  3 

where  At  is  now  (t  -  1950.0)/100  and  the  initial  epoch  depen¬ 
dence  of  the  equatorial  precessional  quantities  (i.e.,  1950.0) 
is  implicitly  included  in  the  above  formulas.  While  P  requires 
six  trigonometric  evaluations  and  many  multiplications  and 
additions  for  its  formation  this  is  overhead  in  the  sense  that 
P  is  formed  once  per  ASTPT  execution  and  not  once  per  minor 
planet.  It  was  the  tabular  source  for  R^  that  ruled  this 
possibility  out  plus  its  relative  inefficiency  compared  to 
#2  or  #3. 

5.  Finally,  since  an  excellent,  simple,  analytical 
20 

form  is  available  for  Rq  the  procedure  used  is  (i)  obtain 

*0  by  Eq.  (8),  (ii)  precess  to  the  epoch  of  date  by  Eqs.  (13) 

to  obtain  n.  -  P4  ,  and  (iii)  then  add  Rrt  which  yields  R. 

—  — o  — (u  — 

C.  The  Location  of  the  Sun 

20 

The  Almanac  for  Computers  provides  Chebyshev  expansions 
for  the  geocentric  right  ascension,  declination,  and  distance 
of  the  Sun  relative  to  the  true  equator  and  equinox  of  date. 
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It  also  provides  similar  series  for  the  nutation  in  longitude 
Aip  and  the  nutation  in  obliquity  Ae .  Accuracy  to  0?6,  3", 

4  x  10  ^  au ,  01'4,  and  0'.'2  respectively  is  available.  A  sub¬ 
routine  (called  SUN)  exists  in  ASTPT  which  computes  these  values 
and  then  corrects  the  geocentric  right  ascension  and  declina¬ 
tion  relative  to  the  true  equator  and  equinox  of  date  (say  Ag, 
Aq)  to  the  geocentric  right  ascension  and  declination  relative 
to  the  mean  equator  and  equinox  of  date  (Ag,  Ag)  by 


A©  38  Ag  -  (cose"  +  sine"sinAgtanAg)  Aip  +  (cosA^tanAg) Ae 
Ag  =  Ag  -  (sine  "cosAg)  Aip  -  (sinAg)Ae 


where  e"  =  e  +  Ae  is  the  true  obliquity  of  date.  From  Rg  = 


R©MAg,  Ag)  the  spherical,  ecliptic,  geocentric  coordinates  for 

« 

the  Sun,  relative  to  the  mean  equator  and  equinox  of  date 
are  computed  by 


cosAgeosAg 

=  |  cosAgSinA^cose  +  sinA^sine 


0 


(14) 


,  -cosA._sinA_sine  +  sinA_cose 
©  ®  Q 


which  is  another  way  of  writing  Eq.  (6).  We  then  obtain  A@  and 

B*  for  use  later. 

© 

All  of  this  is  done  separately  for  the  three  times  and 
represents  overhead. 


D. 


Planetary  Aberration 


R  =  R£/A,  A)  is  the  geometric ,  spherical,  equatorial 

location  for  the  mean  equator  and  equinox  of  date  at  some 

instant  of  interest  (say  12  A.M.).  However,  because  the  speed 

of  light  c  (in  vacuo)  is  not  infinite  and  the  minor  planet  is  a 

reasonable  distance  ('v-  2  A.U.)  away,  the  light  reflected  at 

12  A.M.  won't  reach  us  for  another  ^  10^  seconds  during  which 

time  the  relative  orientation  of  the  Earth  and  the  asteroid 

will  have  changed  (at  something  less  than  l"/min  for  most  main 

belt  asteroids)  by  <  17".  This  is  the  planetary  aberration  and 

we  must  correct  for  it.  Given  the  ^  1'  accuracy  in  the 

position  (cf.  § §VA)  we  can  afford  a  very  simple  form  of 

correction.  In  particular  I  assume  that  the  motion  of  the 

Earth  and  that  of  the  asteroid  is  uniform  and  rectilinear  in 

»  • 

the  time  interval  R/c.  Then,  if  A  and  A  are  the  rates  of 

change  of  A  and  A,  the  values  of  A  and  A  corrected  for  planetary 

•  • 

aberration  are  A  -  AR/c  and  A  -  AR/c.  Alternatively  I  could've 
predated  the  output  time  by  R/c. 

E.  Other  Effects 

Other  effects  might  include  nutation,  annual  aberration 
(both  circular  and  e-terms) ,  diurnal  parallax,  diurnal  aberra¬ 
tion,  astronomical  refraction,  and  parallactic  refraction. 

These  can  easily  exceed  2'  at  a  zenith  distance  of  60°.  How¬ 
ever,  they  are  clearly  not  needed  to  point  the  telescope  (cf. 

§ VA)  and  have  a  negligible  (0.03%)  probability  of  causing  us 
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to  not  provide  output  for  a  minor  planet  at  Z  =  60°01'  which  these 
effects  (principally  astronomical  refraction)  might  have  raised 
up  to  Z  =  59°59'.  They  are  also  not  needed  in  ASTID  because 
the  input  to  ASTID  will  be  the  output  from  the  calibration 

g 

software  PLC  or  SSC(SAO)  at  the  ETS .  These  provide  positions 

relative  to  the  mean  equator  and  equinox  (uncorrected  for 

diurnal  parallax  or  planetary  aberration)  at  the  beginning 

* 

of  the  current  calendar  year  .  Therefore,  if  all  of  these 
small  corrections  were  added  into  the  ASTPT  software,  I'd  just 
have  to  delete  them  to  change  ASTPT  into  ASTID. 

F.  Brightness 

The  most  recent  and  authoritative  summary  of  asteroid 
absolute  B  magnitudes,  phase  functions,  and  the  opposition 
effect  is  in  reference  21.  ASTPT  uses  slightly  older  values 
which  can  result  in  a  difference  of  up  to  O^S  (at  quadrature 
due  to  phase  effects;  at  most  O^IS  near  opposition) .  The 
apparent,  exo-atmospheric  B  magnitude  for  a  minor  planet  r  A.U. 
from  the  Sun,  R  A.U.  from  the  Earth,  and  at  phase  angle  f 
(and  away  from  opposition)  is 

B(rR,f)  =  B  (1 ,  0)  +  51og(rR)  +  Ff  (15) 

where  f  is  the  phase  angle  [see  Fig.  3  and  Eq.  (17)]  and  F  (old) 

=  0I?023/deg.  F(new)  =  01?039/deg.  The  quantity  Ff  is  the 
* 

This  is  a  U.S.  Air  Force  desire. 
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phase  function  in  magnitudes  (and  this  notation  is  not  standard) . 
Equation  (15)  is  valid  for  8°  <  f  <  22°  (quadrature).  When  the 
phase  angle  is  near  zero  (f  =  0  implies  opposition)  the  asteroids 
brighten  perceptibly  faster  than  the  linear  form  in  Eq.  (15) 
would  indicate.  This  is  known  as  the  opposition  effect  and  it 
is  given  by  (old) 

B  (1 ,  f )  =  0*?188|f  |0,545  +  8  °F  -  0*1*584  +  B(1,0) 

(16) 

=  O^ISS { f I  0  -  545  _  0^400  +  B  (1 , 0 ) 

The  new  form  is 

B  (1 ,  f )  =  0*1*134  If  |0*714  +  7°F  (new)  -  0*1*538  +  B(1,0) 

=  0*1*134  |f  |0,714  -  0*1*265  +  B  (1 , 0 ) 

Equation  (16)  is  used  when  |f  |  8°.  Figure  3  implies 

cosf  =  (R2  +  r2  -  R2 ) / (2Rr)  (17) 


Atmospheric  extinction  is  proportional  to  the  air  mass 
involved  and  the  air  mass  is  (roughly)  proportional  to  the 
secant  of  the  (topocentric ,  refracted)  zenith  distance. 

Ignoring  several  small  effects  (astronomical  refraction, 
diurnal  parallax,  nutation,  annual  aberration,  the  deflection 
of  the  vertical,  the  angle  of  the  vertical,  and  color  effects), 
all  of  which  sum  to  <  0*1*1,  we  can  approximate  the  extinction, 


e  =  0?25  secZ 

Z  is  computed  from  Eq.  (18). 

cosZ  =  sinAsimj)  +  cosAcos^cos  (t  -  A)  (18) 

Finally  then,  our  estimate  (good  to  ^  01?!)  of  the 
apparent  S-20  magnitude  of  the  minor  planet  is 

m2Q  =  B(rR,f)  +  0(|f|)  +  e  -  O^S  (19) 

where  0(|f|)  is  the  opposition  effect  evaluated  from  Eq.  (16) 
when  appropriate  and  the  constant  represents  the  S-20/B  color 
term  from  reference  22  using  a  B-V  =  O1?? 7  (see  §11).  The 
scatter  in  B-V  can  introduce  another  +  O^OS  error. 

G.  Visibility 

For  the  purposes  of  ASTPT  an  asteroid  is  "visible"  if  and 
only  if  its  orbital  element  set  (epoch  1950.0,  date  of 
osculation)  passes  the  user  requested  inclination  and  eccen¬ 
tricity  tests,  has  a  value  of  n^p  (evaluated  at  midnight)  from 
Eq.  (19)  that  satisfies  the  user  requested  magnitude  limits, 
and  has  a  value  of  Z  from  Eq.  (18)  less  than  60°  (evaluated 
at  midnight) .  If  an  asteroid  passes  all  of  these  tests  then 
the  6  P.M.  and  6  A.M.  positions  are  computed,  the  geocentric 
angular  velocity  computed  at  midnight  from  the  three  positions 
by  differentiating  a  three-point  Lagrange  interpolation  formula, 
the  planetary  aberration  correction  applied,  the  time  to  move 
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7 '.'5  calculated,  output  prepared,  and  the  requisite  information 
stored  for  a  potential  sort. 

To  paraphrase  the  Explanatory  Supplement  to  The  Astronomical 
Ephemeris  and  the  American  Ephemeris  And  Nautical  Almanac,  even 
when  the  computations  are  not  performed  as  discussed  above  or 
as  in  § VI ,  neither  are  they  done  very  differently  nor  is  any 
appreciable  (at  +1')  error  introduced. 
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VI .  THE  TRICKS 


Clearly  the  computation  of  a  spherical,  equatorial, 
geocentric  position  of  an  asteroid  for  an  arbitrary  mean 
equator  and  equinox  represents  a  fair  amount  of  calculation. 

To  make  ASTPT  and  especially  ASTID  run  as  quickly  as  possible 
this  should  be  avoided  whenever  feasible.  My  methods  of 
doing  this  are  given  in  subsection  D  below.  A  very  simple 
way  of  saving  CPU  time  is  to  run  in  single  precision  (21 
bits  =  0  71  on  my  computer)  instead  of  double  precision 

4r 

(37  bits  -  1  y  arc  sec) .  All  overhead  is  performed  in 
double  precision  and  then  truncated.  All  other  computations 
are  performed  in  single  precision. 

A.  Heliocentric  Locations 

The  vector  for  midnight  is  calculated  as  outlined  in 
§VA  (i.e.,  rigorously).  The  6  P.M.  and  6  A.M.  values  are  then 
computed  differentially  via 

A  (t  +  At)  =  K  (t)  +  (dA  (t) /dt ) At  (20a) 

— o  —o  — o 

where  t  =  midnight.  At  =  +6*\  and,  from  Eqs.  (1,  8) 
dA  (t) 

—  =  n(  (cosEJE^  -  (sinEjA^J/U  -  ecosE)  (20b) 

* - 

Overhead  consists  of  the  creation  of  the  master  asteroid  file 
(e.g.,  A  ,  B,  n,  etc.),  the  geocentric  location  of  the  Sun,  the 
local  mean  sidereal  time,  the  general  precession  matrix  P, 
and  the  ecliptic  latitude  and  longitude  of  the  astronomical 
zenith. 
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where  it's  the  12  A.M.  value  of  the  eccentric  anomaly  that  is 
used.  The  second-order  error  can  scarcely  exceed  0'.'5  for  an 
object  moving  at  l"/min. 

B.  General  Precession  and  Parallax 

The  vector  >i  for  midnight  is  rigorously  computed  from  Eqs. 

(13)  (i.e.,  fi  =  P^  with  the  values  for  the  equatorial  preces- 

sional  quantities  evaluated  at  midnight).  From  Eq.  (8)  we  see 

that  this  involves  the  computation  of  PA  and  PB  .  If  we  don't 

— o  — o 

h 

adjust  P  for  the  +  6n  and  use  Eqs.  (20),  then  no  additional 
matrix  multiplications  are  required  to  obtain  the  6  A.M.  or  6  P.M. 
values  of  n.  relative  to  the  mean  equator  and  equinox  of  midnight. 
Since  m  -  46'.'l/yr  and  n  =  20'.'/yr,  the  maximum  error,  from  Eqs. 

(11,  12)  is  <  0'.'05  for  any  place  on  the  celestial  sphere.  And  this 
is  the  greatest  of  all  upper  bounds. 

As  the  location  of  the  Sun  is  highly  accurate  and  Fig.  3 
is  rigorous,  the  parallax  (heliocentric-geocentric)  correction 
is  performed  rigorously  at  the  three  times  of  interest.  Moreover, 
rigorous  parallax  for  ecliptic  coordinates  can  be  written  as 


tan  (A  -  X)  = 


cosBrtsin(Art  -  X) 

_ U _ <y _ 

(r/R0)cosB  +  cosBqcos(Aq 


(21a) 
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where 


tany  =  tanBQcos[(A  -  A)/2]sec[AQ  -  (A  +  A)/2]  (21c) 

This  will  be  used  later  (§VID) .  Note  though,  that  since 
|Bq|  <_  1"2  unless  A@  -  (A  +  A)/2  =  +  90°,  y  =  0  too. 

I 

C.  Planetary  Aberration 

The  only  approximation  made  here,  not  already  detailed 
in  §VD,  is  to  use  the  same  planetary  aberration  correction  at 
all  three  times  of  interest.  This  error  is  comparable  to  that 
made  in  the  rectangular,  heliocentric,  1950.0  locations  because 
I  it  depends  on  the  relative  acceleration  of  the  Earth  and  the 

asteroid  over  the  +  6h. 

D.  Visibility 

An  asteroid  must  pass  four  levels  of  tests  before  its 
position  is  ever  computed.  Two  more  levels  of  tests  are 
required  after  the  position  is  available.  These  concern  zenith 
distance  and  apparent  magnitude  (§VID6).  Additional,  preliminary 
zenith  distance  (§VID4)  and  apparent  magnitude  (§VID2)  tests 
are  also  performed. 

1.  User  Inclination  and  Eccentricity  Limits 
Obviously  the  quickest  methods  of  elimination  are  to  ensure 

that  i^  £  iQ  £  iu  and  e^  <  e  <_  eu  and  these  are  performed  first. 

2.  User  Apparent  Magnitude  Test  (Preliminary) 

For  the  huge  majority  of  the  asteroids  in  the  main  belt 

the  perihelion  opposition  magnitude  [m^>;  r  =  a(l  -  e)  =  q, 
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R  =  q  -  1J  is  brighter  than  the  mean  opposition  magnitude 
(r=a,  R=a-1)  which  is  brighter  than  the  aphelion 
opposition  magnitude  [r  =  a (1  +  e)  ,  R  =  a (1  +  e)  -  1]  which, 
in  turn,  is  brighter  than  the  mean  quadrature  magnitude 
[m  ;  r  =  a,  R  =  (a2  -  l)^2,  f  =  sin  ^(1/a)].  Including 
extinction  for  one  air  mass  these  are  (in  the  blue) 


m^  =  B  (1 , 0 )  +  0  (o)  +  51og[q(q  -  1)]  +  0?25 

mn  =  B  (1 , 0 )  +  51og [a  (a2  -  1)1/2]  +  Fsin-1(l/a)  +  0?25 


(22) 


Therefore,  if  m^  <  m^  or  m^,  >  m^,  then  the  minor  planet  can 
be  eliminated  (when  itu  and  m  are  corrected  for  the  S-20/B 

O  I  i 

color  term  as  in  Eq.  (19)].  The  value  of  the  perihelion 
opposition  magnitude  and  mean  quadrature  magnitudes,  so 
adjusted,  are  on  the  master  asteroid  tape  file  and  read  in 
with  the  orbital  element  set,  A^,  8q,  etc. 

3.  Right  Hemisphere  Test 

* 

Roughly  half  of  all  of  the  asteroids  are  going  to  be  in 
the  wrong  hemisphere  of  the  celestial  sphere  at  any  given  time. 
If  they  can  be  efficiently  eliminated  from  consideration  early 
then  a  great  deal  of  computation  will  be  saved.  This  test  is 
designed  to  do  just  this.  Consider  Fig.  4.  Only  those  minor 


* 

This  assumes,  among  other  things,  a  uniform  distribution  of 
mean  longitudes.  This  is  not  the  case,  primarily  due  to  per¬ 
turbations  from  Jupiter.  There  are  noticeable  peaks  in  the 
longitude  of  ascending  node  distribution  function  near  H-vj  and 
ftV|  +  180°.  A  similar,  smaller  effect,  is  present  in  the 
argument  of  perihelion. 
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4.  Minor  planets  with  geocentric  ecliptical  longitudes. 


f 


planets  with  geocentric  ecliptic  longitudes  Ae (Ag  +  90°,  A^ 

+  270°)  are  of  interest.  But  it  is  Xq  we  can  get  from  the 
orbital  element  set,  not  A  or  even  X. 

We  start  from 

X  =  ft  +  tan  ^  (cosi  tanu  )  (23) 

o  o  o  o 

where  uq  is  the  true  argument  of  latitude  =  u)q  +  v.  (v  is  the 

* 

true  anomaly  we  didn't  compute  in  §VA. )  Since  ,  by  and  large, 

e  <0.3  and  i  <  30°  we  can  use 
o  —  o  — 

Ao  ~  Ao  =  fio  +  uo  +  M  +  2esinM  +  (5e2/4)sin2M  (24) 

where  M  is  the  mean  anomaly,  the  Taylor  series  of  Eq.  (23)  about 
iQ  =  0  has  been  used  through  first-order  terms  in  i,  and  the 
Taylor  series  of  Eqs.  (1,3)  about  e  =  0  have  been  used  through 
second-order  terms  in  e.  Numerically*,  the  worst  Eq.  (24)  is  is 
2?7  +  1?6  averaged  over  a>Q  =  0(15)360°,  M  =  0  (15)360°  for  test 
values  of  iQ  =  0(5)30°,  e  =  0(0.05)0.30.  Not  surprisingly  this 
occurs  for  iQ  =  30°  and  e  =  0.30.  Dropping  the  quadratic  term 
in  e  raises  2?7  +  1?6  to  9f2  +28?6.  Now  that  we’ve  got  a  +5° 

_ 

Compare  with  reference  17. 

+ 

These  results  refer  to  the  mean  value  of  | X  -  X ' |  and  its 
standard  deviation  about  the  mean,  i.e. ,  mea8  | X  °—  X  * |  at  a 

u)  =  0(15)360°  ° 

0  =  0  (15)360° 

fixed  i  and  e  is  computed  and  2?7  was  the  largest  value  in  the 
iQ,  e  ranges. 
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approximation  for  XQ  we  need  to  precess  it  to  the  equator 
and  equinox  of  date.  In  terms  of  the  quantities  a,  b,  and  c 
introduced  in  Eq.  (9)  general  precession  for  spherical 
ecliptic  coordinates  is  given  by  (to  the  second-order) 

\  =  X  +  a  -  bcos(X^  +  c)tan3 
o  o 

3  =  3Q  +  bsin(Xo  +  c) 

But  a  =  0?46 ,  b  a  15 ?5  and  since  3  is  small  (by  hypothesis) 
we  can  ignore  general  precession  completely  in  both  X  and 
3  (at  the  +5°  level).  Just  to  be  safe  the  general  precession 
in  longitude  has  been  kept  in  (it's  overhead)  so 

X  =  X '  =  X '  +  a  a  X  +  a 
o  o 

To  get  A  (or  an  approximation  to  it)  we  look  at  Eq.  (21a) 
and  realize  that  BQ  *  0  (within  172).  Note  also  that  R0 
=  1  A.U.  (within  e0  =  0.01675).  Therefore, 

sin  (A_  -  X) 

tan  (A  -  X)  =  - 

(r/R0)cos3  +  cos (A0  -  X) 

To  at  worst  1?2  +  1?5  [for  (r/R0)cos3  *  1.5(0. 5)5.0  A.U. 
and  averaged  over  A&  -  X  =  0(15)360°  as  above,  e.g.,  for  the 
mean  of  | A  —  X | ]  we  can  approximate  this  by 
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A  -  X  =  A 1 


sin (Aq  -  A ' ) 


-  A'  = 


(r/RQ)cosB  +  cos(Aq  -  A') 


sin (Aq  -  A*) 

(r/R@)cosB  +  cos  (A@  -  A') 


Finally,  for  rcosg  we  use  a (1  -  ecosM)  which  is  accurate  to 
0.066  +  0.047  in  (r/a)cosB  (averaged  over  u)Q,  M;  the  worst  here 
is  for  i  =  30°,  e  =  0) . 

We  now  have  an  approximate  (<  5°)  value  for  A  at  the 
cost  of  sinM,  cosM,  sin(A0  -  A'),  and  cos (AQ  -  A').  And  none 
of  these  are  wasted  (cf.  §§VID4,  VID5).  Using  our  approximate 
value  for  A  we  reject  all  minor  planets  not  within  95°  of 
(AQ  +  180° )mod (360° ) .  If  either  e  >  0.3  or  i  >  30°  then  this 
computation  is  bypassed. 

4.  Above  Horizon  Test  (Preliminary) 

Only  one-quarter  of  the  celestial  sphere  is  above 
altitude  =  30°.  Hence,  if  we  can  eliminate  those  minor 
planets  below  this  (and  do  it  cheaply) ,  another  large  savings 
will  be  made.  Of  course,  it  won't  be  three-quarters  because 
half  of  the  celestial  sphere  has  already  been  ruled  out.  I 
expect  that  about  half  of  the  remaining  candidates  will  be 
eliminated  here.  To  do  this  though  we  need  to  compute  the 
minor  planet's  zenith  distance.  To  use  Eq.  (18)  for  cosZ  clearly 
involves  knowledge  of  the  equatorial  position.  Hence  we  need 
another  method.  This  trick  consists  of  computing  cosZ  in  the 
ecliptic  coordinate  system. 


Jr 


If  B  ,  L  are  the  spherical,  geocentric,  mean  equator 
z  z 

and  equinox  of  date,  ecliptic  latitude  and  longitude  of  the 
astronomical  zenith  and  B,  A  the  corresponding  coordinates 
for  the  asteroid,  then 


cosZ  =  £  (A , B)  .  £ (L  , B  ) 

z  z 


(26) 


Since  Eq.  (14)  represents  the  general  relationship  between 
spherical  ecliptic  and  spherical  equatorial  coordinates, 
and  the  declination  of  the  astronomical  zenith  is  <{>  while 
its  right  ascension  is  just  t. 


cos<J>cost 

£(L  ,B  )  =  I  cos<J>sinTcose  +  sin<)>sine 


-cos<|>sinTsin£  +  sin<|>cose 


(27) 


/ 


Note  that  this  is  exact  and  £ (L  , B  )  represents  overhead. 

z  z 

The  approximation  for  B  is  straightforward.  We  have, 
rigorously. 


B_  =  sin-1  (sini  sinu  ) 
o  o  o 


(28) 


whence 


B0  ~  Pq  =  i0sin(o>0  +  M  +  2esinM)  (29) 

never  worse  than  1?3  +  1?0  [cf.  the  discussion  following 

Eq.  (24);  again  for  iQ  30°,  e  0.3  and  this  is  the  30°,  0.3 

value].  We  already  know  that  general  precession  is  negligible, 


w 


so,  with  this  accuracy  B  =  B^.  Parallax  is  also  simple. 
From  Eqs .  (21),  with  sinB0cscy  -  B0/y  -  tanB0/tany 

=  cos  [  A0  -  (A  +  A)/2]  sec  [  (A  -  A)/2], 


tan(B  -  B) 


-sinB _ 

(y/BQ) (r/R0)  +  cosB 


because  y  =  0  unless  A0  -  (A  +  A)/2  =  +  90°.  And  if  we 
neglect  parallax  for  the  moment,  this  is  the  same  as 
AQ  -  A  =  +  90°  so  we  realize  that  the  points  of  greatest 
trouble  have  already  been  declared  uninteresting.  Finally, 
then 


-  sinB' 

B-B  -  B '  -8*  -  - 

(y/B0) {r/R©)  +  c°s&' 


_ -  sinB' _ 

(y/B0) (r/R0)  +  cosB ' 


(30) 


with  an  accuracy  of  at  worst  2 ? 2  +  1?4  [averaged  over  B' 

=  -90(10)90°  for  (y/B0)  (r/R0)  =  -5  (0 . 5) -1. 5] .  Substituting 

2 

for  B'  using  Eq.  (29)  and  using  r/a  =  1  -  ecosM  +  (e  /2). 

(1  -  cos2M)  (good  to  0.0003  +  0.0019  even  at  e  =  0.3)  completes 
the  computation  of  B' .  The  secant  of  (A'  -  A')/2  is  computed 
from  1  +  (A'  -  A')2/8  since  |A  -  A|  <  20°.  The  quantity 
* 

From  Eq.  (21a)  wjth  BQ  =  0  R  =  1  A.U.  it  follows  that 
| tan (A  -  A) |  <.  cos  B  -  1)  '  (this  is  rigorous).  At  rcosB 

=  2.7  A.U.  this  implies  jA  -  A|  <  22°. 
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cos [Aq  -  (X*  +  A*)/2J  =  cos [Aq  -  X'  +  (X'  -  A')/2]  is 
calculated  as 

cos  (A0  -  X  * )  [1  -  (X'  -  A '  )  2/8]  -  (-X  '  ~  A  - }  sin(A0  -  X') 

Therefore  only  the  trigonometric  function  in  Eq.  (29)  is 
needed . 

We  started  by  computing  four  trigonometric  functions  in 
the  right  hemisphere  test.  This  is  less  work  than  just  solving 
Kepler's  equation  (see  below).  Hence,  if  the  minor  planet 
can  be  eliminated  on  this  basis  a  huge  savings  will  result. 

Of  those  four  trigonometric  functions  two  are  used  in  this 
test.  And  if  the  minor  planet  is  eliminated  from  further 
consideration  now,  then  a  total  of  five  trigonometric 
functions  have  been  computed.  This  is  still  cheaper  than 
solving  Kepler's  equation.  The  actual  test  used  is  that  if 
cosZ,  computed  from  Eq.  (26)  using  Eqs.  (25,  30)  for  A  and  B 
is  <  1/3  (which  corresponds  to  the  error  in  cosZ  from  the 
accumulation  of  all  of  the  various  approximates,  i.e.,  at 
worst  3?9  +  2?1  in  A  and  3?5  +  1?8  in  B)  then  the  minor 
planet  is  rejected. 

5.  Kepler's  Equation 

At  this  point  I'm  out  of  tricks  and  the  full  midnight 
position  must  be  computed.  As  discussed  in  §VA  the  first 
step  is  to  compute  the  mean  anomaly  M  (already  done)  and 
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then  solve  Kepler’s  equation,  Eq.  (1).  Four  such  methods  have 
been  tried  over  the  e  =  0(0.05)0.30,  M  =  0(15)360°  ranges.  The 
fastest  is  Newton's  method. 


E  =  E  .  , 
new  old 


(Eold  ~  esinEold  "  M) 


1  -  ecosE 


old 


if 


new 


-  E  , , I  is  small  enough,  then  E  =  E  and  terminate 
old 1  3  new 


otherwise  E  ,  .  =  E  and  repeat 
old  new 


with  the  original  Eo^  =  M  +  esinM.  On  the  average  this  requires 
2.5  iterations.  The  second  fastest  is  Newton's  method  with  a 
starting  guess  of  Eq^  =  M  (3  iterations  on  average)  .  But  for 
this  case  the  sine  and  cosine  of  M  are  needed,  we  already  have 
these,  and  Newton's  method  starting  with  E^^  =  M  and  the  addi¬ 
tional  headstart  of  sinM,  cosM  is  the  fastest  of  all.  Hence, 
when  passing  through  the  right  hemisphere  and  zenith  distance 
tests  only  a  single  trigonometric  function,  the  one  in  Eq.  (29), 
is  "wasted". 

6 .  The  Rest 

Any  asteroid  getting  to  the  solution  of  Kepler's 
equation  (and  this  automatically  includes  the  few  for  which 
either  e  >  0.3  or  iQ  >  30°)  has  its  spherical,  equatorial, 
geocentric ,  mean  equator  and  equinox  of  date  location  computed 
as  outlined  in  §§VA,  VB,  VC,  VE,  VIA,  and  VIB.  Then  cosZ  is 
calculated  from  Eq.  (18).  If  it's  <_  1/2  the  minor  planet 
is  rejected.  If  not,  then  it^q  is  calculated  from  Eq.  (19) 


and  tested  against  mf.  If  it  passes  this  test,  then 
the  6  P.M.  and  6  A.M.  locations  are  computed,  the  geocentric 
angular  speed  calculated  {A(12  A.M.)  =  2[A(6  A.M.)  -  A(6  P.M.)], 
A  per  hour,  etc.},  the  planetary  aberration  applied,  the  time 
to  move  7 '.'5  computed,  the  position  angle  and  direction  of  the 
geocentric  angular  velocity  obtained,  output  readied,  and 
much  of  this  stored  if  a  sort  has  been  requested. 

E.  ASTID  Visibility 

In  this  mode  only  asteroids  within  a  few  degrees  of  a 
particular  geocentric  ecliptic  latitude  and  longitude  are 
of  interest.  From  the  equatorial  coordinates  of  the  asteroid 
being  observed  we  compute  the  corresponding  ecliptic 
coordinates  [by  inverting  Eq.  (6)  with  appropriately  accurate 
general  precession  and  planetary  aberration  included;  we 
started  with  a  geocentric  right  ascension  and  declination 
for  the  mean  equatorial  equinox  of  1980.0  (for  1980)].  Call 
these  L* ,  B* .  For  each  asteroid  on  file  we  calculate  A'  from 
Eq.  (25)  and  reject  it  as  a  candidate  for  identification  if 
A'  -  L*  \  ~>  5°.  On  the  average  this  eliminates  35/36  =  97.2% 
of  the  entire  catalog.  If  a  minor  planet  passes  this  step 
we  get  B'  as  in  Eq.  (30)  and  then  eliminate  it  if  | B '  -  B* j  > 

5°.  This  removes  170/180  of  the  remainder,  leaving  about  3 
(plus  the  high  e,  i  )  for  further  consideration.  Of  course, 
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VII.  ASTID 


The  ASTID  visibility  tricks  given  in  §VIE  are  a  natural 
outgrowth  of  those  used  in  ASTPT.  We  can  obtain  a  rough 
estimate  of  their  relative  speeds  of  execution  as  follows: 
Assume  that  the  eccentricity  and  inclination  limits  are  at 
their  default  values.  If  the  right  hemisphere  test  eliminates 
half  of  these  and  the  preliminary  above  the  horizon  test 
eliminates  half  the  remainder,  then  we  have  ^  550  left.  In 
addition  we  have  'v  50  high  a  or  high  i  objects  so  ASTPT  must 
fully  process  ^  600  minor  planets.  In  the  ASTID  outlined 
above  we  have  ^  (1  -  35/36)  (1  -  17/18)2200  +  50  to  process, 

or  53.  Thus,  assuming  that  the  program  overhead  is  negligible 
compared  to  the  location  computation,  the  relative  speed  should 
be  ^  12.  A  bit  more  conservative  estimate  is  the  ^  10  quoted 
in  the  Introduction.  This  is  still  far  too  long  6  seconds) . 

In  particular,  when  ASTPT  is  run  during  the  day,  for 

•  • 

that  evening,  we  will  save  the  midnight  A,  A,  A,  and  A  values 
calculated  for  all  asteroids  that  were  visible.  Now,  from 
the  current  telescope  position  and  time  (i.e.,  right  ascension 
and  declination)  ASTID  computes  the  corresponding  asteroid 
position  direction  cosines  using  A  =  A (12  A.M.)  +  AAt,  etc. 
and  we  investigate  the  size  of 

cosD  =  i(A,A)  •  MAt,At) 

where  A, A  are  for  a  particular  asteroid  at  this  instant  and 


At,At  are  for  the  telescope  position.  D  is  the  angular 

-20 

distance  between  the  two  points  and  we  test  for  cosD  <1-2 
or  D  <  4!7.  This  includes  orbital  element  set  errors,  telescope 
errors,  and  computational  errors  (there  are  no  errors  in 
A, A  due  to  the  parallax,  general  precession,  or  other  approxi¬ 
mations  used  elsewhere  in  ASTPT) .  Clearly  ASTID  is  now 

•  • 

input-output  bound  (i.e.,  reading  the  A, A, A, A  file).  The 
addition  of  this  step  to  ASTPT  will  slow  it  down  negligibly. 


1 


APPENDIX:  User's  Guide 

ASTID  is  a  real  time  program  at  the  ETS  and  it  is  executed 
by  pushing  one  of  the  buttons  on  the  operating  console  (i.e., 
the  one  labeled  ASTID) .  The  execution  of  ASTPT  is  started 
by  performing  a  console  interrupt  and  then  typing  ASTPT/ /ACT, ,TM. 
Sample  sessions  of  the  user  interaction  with  ASTPT  are  shown 
in  Tables  A1  -  A5 . 

Table  A1  illustrates  the  full,  normal  method  of  choosing 
the  default  options.  In  this,  and  the  other  four  tables,  user 
inputs  are  underlined.  All  dialog  appears  on  the  line  printer 
and  on  the  CRT.  When  the  code  executes  a  top  of  form  command 
it's  been  indicated  by  a  dashed  line  here  (to  save  space) . 

Table  A2  illustrates  the  full,  normal  method  of  entry 
with  several  quantities  being  changed.  The  notation  is  fairly 
obvious  (m^  ->  MB,  n^g  -*■  MASTEROID,  iQ-*-  INCL,  a  -*■  RA,  etc.) 

Table  A3  illustrates  the  form  and  nature  of  the  error 
messages. 

Table  A4  illustrates  an  attempt  at  the  quick  method  of 
entry.  This  failed  and  the  user  chose  to  return  to  the  full, 
normal  mode.  To  save  space  this  is  not  given  here. 

Table  A5  illustrates  a  successful  quick  mode  session. 
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TABLE  A1 


ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 

IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 

2 _ 5 

YOU  ENTERED  2  5 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO  *  IF  NOT  THEN 
ANY  INTEGER 
-9 


OUTPUT  IS  FOR  2  5  80  AT  6PM,  12AM,  AND  6AM  LOCAL 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 
55 

DO  YOU  WISH  TO  IMPOSE  MAGNITUDE  LIMITS? 

THE  DEFAULTS  ARE  0  =  MBRIGHT  .LE.  MASTEROID  .LE.  MFAINT  =  18.5 

TYPE  A  ONE  TO  ONLY  CHANGE  MBRIGHT 

TYPE  A  TWO  TO  ONLY  CHANGE  MFAINT 

TYPE  A  THREE  TO  CHANGE  BOTH  MBRIGHT  AND  MFAINT 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MBRIGHT  NOR  MFAINT 

4 

YOU  ENTERED  4 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN  ANY 
INTEGER 
4 
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THE  MAGNITUDE  LIMITS  ARE  0.00  .LE.  MASTERQID  .LE.  18.50 

DO  YOU  WISH  TO  IMPOSE  INCLINATION  LIMITS? 

THE  DEFAULTS  ARE  0  =  MIN  INCL  .LE.  INCL  .LE  MAX  INCL  =  180 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  INCL 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  INCL 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  INCL  AND  MAX  INCL 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  INCL  NOR  MAX  INCL 

4 

YOU  ENTERED  4 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN  ANY 
INTEGER 
4 


THE  INCLINATION  LIMITS  ARE  0.00  .LE.  INCL  .LE.  180.00 

DO  YOU  WISH  TO  IMPOSE  ECCENTRICITY  LIMITS? 

THE  DEFAULTS  ARE  0  =  MIN  ECC  .LE.  ECC  .LE  MAX  ECC  =  l 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  ECC 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  ECC 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  ECC  AND  MAX  ECC 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  ECC  NOR  MAX  ECC 

4 

YOU  ENTERED  4 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
4 
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THE  ECCENTRICITY  LIMITS  ARE  0.00  .LE.  ECC  .LE.  1.00 

ENTER  TAPE  UNIT  »  --  1  FOR  MTl  OR  2  FOR  MT2 
1 

YOU  ENTERED  1 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN  ANY 
INTEGER 

1 


TAPE  SHOULD  BE  MOUNTED  ON  MTl 

THE  OUTPUT  IS  AUTOMATICALLY  GIVEN  BY  INCREASING  ASTEROID  NUMBER 

YOU  MAY  HAVE  IT  SORTED  BY  INCREASING  RA,  MAG.  OR  ANGULAR  SPEED 

THESE  SORTS  REQUIRE  SOME  EXECUTION  TIME  AND  YOU  HAVE  EIGHT  CHOICES 

TYPE  A  ZERO  FOR  NO  ADDITIONAL  SORTING 

TYPE  A  ONE  FOR  AN  RA  SORT 

TYPE  A  TWO  FOR  A  MAG  SORT 

TYPE  A  THREE  FOR  AN  ANGULAR  SPEED  SORT 

TYPE  A  FOUR  FOR  AN  RA  AND  A  MAG  SORT 

TYPE  A  FIVE  FOR  AN  RA  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SIX  FOR  A  MAG  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SEVEN  FOR  ALL  THREE  SORTS 


YOU  ENTERED  0 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN  ANY 


INTEGER 


YOU  DESIRE  NO  ADDITIONAL  SORTING 

EXECUTION  BEGINS  N 
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TABLE  A2 


ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 
IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 
5  2 

YOU  ENTERED  5  2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

88 


OUTPUT  IS  FOR  5  2  80  AT  6PM,  12AM,  AND  6AM  LOCAL 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 

777 

DO  YOU  WISH  TO  IMPOSE  MAGNITUDE  LIMITS? 

THE  DEFAULTS  ARE  0  =  MBRIGHT  .LE.  MASTEROID  .LE.  MFAINT  =  18.5 

TYPE  A  ONE  TO  ONLY  CHANGE  MBRIGHT 

TYPE  A  TWO  TO  ONLY  CHANGE  MFAINT 

TYPE  A  THREE  TO  CHANGE  BOTH  MBRIGHT  AND  MFAINT 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MBRIGHT  NOR  MFAINT 

1_ 

YOU  ENTERED  1 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

1_ 

ENTER  DESIRED  MBRIGHT 
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12.56 


YOU  ENTERED  12.56 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO  »  IF  NOT  THEN 
ANY  INTEGER 
3 


THE  MAGNITUDE  LIMITS  ARE  12.56  .LE.  MASTERO ID  .LE.  18.50 

DO  YOU  WISH  TO  IMPOSE  INCLINATION  LIMITS? 

THE  DEFAULTS  ARE  0  =  MIN  INCL  .LE.  INCL  .LE  MAX  INCL  =  180 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  INCL 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  INCL 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  INCL  AND  MAX  INCL 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  INCL  NOR  MAX  INCL 

2. 

YOU  ENTERED  2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 

J5 

ENTER  DESIRED  MAXIMUM  INCLINATION  IN  DEGREES 
23 

YOU  ENTERED  23.00 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 
23 
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THE  INCLINATION  LIMITS  ARE  0.00  .LE.  INCL  .LE.  23.00 
DO  YOU  WISH  TO  IMPOSE  ECCENTRICITY  LIMITS? 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  ECC 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  ECC 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  ECC  AND  MAX  ECC 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  ECC  NOR  MAX  ECC 

3 

YOU  ENTERED  3 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN  ANY 
INTEGER 


3 

ENTER  DESIRED  MINIMUM  ECCENTRICITY 

0.01 

YOU  ENTERED  0.01 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
3 

ENTER  DESIRED  MAXIMUM  ECCENTRICITY 
.99 

YOU  ENTERED  0.99 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 


ANY  INTEGER 


THE  ECCENTRICITY  LIMITS  ARE  O.Oi  .LE.  ECC  .LE.  0.99 
ENTER  TAPE  UNIT  0  —  1  FOR  MTj  OR  2  FOR  MT2 
2 

YOU  ENTERED  2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

2 


TAPE  SHOULD  BE  MOUNTED  ON  MT 2 

THE  OUTPUT  IS  AUTOMATICALLY  GIVEN  BY  INCREASING  ASTEROID  NUMBER 
YOU  MAY  HAVE  IT  SORTED  BY  INCREASING  RA,  MAG,  OR  ANGULAR  SPEED 
THESE  SORTS  REQUIRE  SOME  EXECUTION  TIME  AND  YOU  HAVE  EIGHT  CHOICES 


TYPE  A  ZERO  FOR  NO  ADDITIONAL  SORTING 

TYPE  A  ONE  FOR  AN  RA  SORT 

TYPE  A  TWO  FOR  A  MAG  SORT 

TYPE  A  THREE  FOR  AN  ANGULAR  SPEED  SORT 

TYPE  A  FOUR  FOR  AN  RA  AND  A  MAG  SORT 

TYPE  A  FIVE  FOR  AN  RA  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SIX  FOR  A  MAG  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SEVEN  FOR  ALL  THREE  SORTS 


YOU  ENTERED  7 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 


7 


YOU  DESIRE  ALL  THREE  -  THIS  WILL  TAKE  AWHILE 


TABLE  A3 


ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 
IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 
13  2 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 
IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 

1  32 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 
IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 
1  1 

YOU  ENTERED  1  1 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

0 

ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 
IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 

2  5 

YOU  ENTERED  2  5 
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IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 


ANY  INTEGER 

666 


OUTPUT  IS  FOR  2  5  80  AT  6PM,  12AM,  AND  6AM  LOCAL 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 
DD 

TRY  AGAIN  -  BAD  READ 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 

6 

DO  YOU  WISH  TO  IMPOSE  MAGNITUDE  LIMITS? 

THE  DEFAULTS  ARE  0  =  MBRIGHT  .LE.  MASTEROID  .LE.  MFAINT  =  18.5 

TYPE  A  ONE  TO  ONLY  CHANGE  MBRIGHT 

TYPE  A  TWO  TO  ONLY  CHANGE  MFAINT 

TYPE  A  THREE  TO  CHANGE  BOTH  MBRIGHT  AND  MFAINT 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MBRIGHT  NOR  MFAINT 

3 

YOU  ENTERED  3 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 
3 

ENTER  DESIRED  MBRIGHT 
19 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  DESIRED  MBRIGHT 
15. 
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YOU  ENTERED  15.00 


IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 

2 

ENTER  DESIRED  MFAINT 
14.99 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  DESIRED  MFAINT 
17.5 

YOU  ENTERED  17.50 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 
3 


THE  MAGNITUDE  LIMITS  ARE  15.00  .LE.  MASTERO ID  .LE.  17.50 

DO  YOU  WISH  TO  IMPOSE  INCLINATION  LIMITS? 

THE  DEFAULTS  ARE  0  =  MIN  INCL  .LE.  INCL  .LE  MAX  INCL  =  180 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  INCL 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  INCL 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  INCL  AND  MAX  INCL 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  INCL  NOR  MAX  INCL 

3 

YOU  ENTERED  3 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 


ANY  INTEGER 


ENTER  DESIRED  MINIMUM  INCLINATION  IN  DEGREES 


Z± 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  DESIRED  MINIMUM  INCLINATION  IN  DEGREES 
5 

YOU  ENTERED  5.00 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
5 

ENTER  DESIRED  MAXIMUM  INCLINATION  IN  DEGREES 
15.67 

YOU  ENTERED  15.67 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 
7 


THE  INCLINATION  LIMITS  ARE  5.00  .LE.  INCL  .LE.  15.67 

DO  YOU  WISH  TO  IMPOSE  ECCENTRICITY  LIMITS? 

THE  DEFAULTS  ARE  0  =  MIN  ECC  .LE.  ECC  .LE  MAX  ECC  ■  1 

TYPE  A  ONE  TO  ONLY  CHANGE  MIN  ECC 

TYPE  A  TWO  TO  ONLY  CHANGE  MAX  ECC 

TYPE  A  THREE  TO  CHANGE  BOTH  MIN  ECC  AND  MAX  ECC 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MIN  ECC  NOR  MAX  ECC 

3 

YOU  ENTERED  3 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO,  IF  NOT  THEN 
ANY  INTEGER 
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ENTER  DESIRED  MINIMUM  ECCENTRICITY 
1.05 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  DESIRED  MINIMUM  ECCENTRICITY 
._5 

YOU  ENTERED  0.50 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

8 

ENTER  DESIRED  MAXIMUM  ECCENTRICITY 
•2 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  DESIRED  MAXIMUM  ECCENTRICITY 
0.51 

YOU  ENTERED  0.51 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
5 

THE  ECCENTRICITY  LIMITS  ARE  0.50  .LE.  ECC  .LE.  0.51 
ENTER  TAPE  UNIT  *  —  1  FOR  MTl  OR  2  FOR  MT2 
3 

AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  TAPE  UNIT  0  --  1  FOR  MTl  DR  2  FOR  MT2 
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AT  LEAST  ONE  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
ENTER  TAPE  UNIT  0  --  1  FOR  MTl  OR  2  FOR  MT2 
2 

YOU  ENTERED  2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

2 


TAPE  SHOULD  BE  MOUNTED  ON  MT2 

THE  OUTPUT  IS  AUTOMATICALLY  GIVEN  BY  INCREASING  ASTEROID  NUMBER 

YOU  MAY  HAVE  IT  SORTED  BY  INCREASING  RA,  MAG.  OR  ANGULAR  SPEED 

THESE  SORTS  REQUIRE  SOME  EXECUTION  TIME.  AND  YOU  HAVE  EIGHT  CHOICES 

TYPE  A  ZERO  FOR  NO  ADDITIONAL  SORTING 

TYPE  A  ONE  FOR  AN  RA  SORT 

TYPE  A  TWO  FOR  A  MAG  SORT 

TYPE  A  THREE  FOR  AN  ANGULAR  SPEED  SORT 

TYPE  A  FOUR  FOR  AN  RA  AND  A  MAG  SORT 

TYPE  A  FIVE  FOR  AN  RA  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SIX  FOR  A  MAG  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SEVEN  FOR  ALL  THREE  SORTS 


YOU  ENTERED  5 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 


5 


YOU  DESIRE  A  RIGHT  ASCENSION  AND  AN  ANGULAR  VELOCITY  SORT 
EXECUTION  BEGINS  NOW 
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TABLE  A4 


ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (PROGRAM  ASTPT  3/80) 

AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 

IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 

2 _ 5 

YOU  ENTERED  2  5 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

2 

OUTPUT  IS  FOR  2  5  80  AT  6PM,  12AM,  AND  6AM  LOCAL 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 

0 

ENTER  MB, MF . IL , IU , EL , EU .TAPE  UNIT  * 

15.5  18  0  14  0  .23  99 
YOU  GOOFED  IT  UP 

AT  LEAST  ON  INPUT  OUT  OF  BOUNDS  -  BE  REASONABLE 
IF  YOU  *  D  RATHER  DO  IT  THE  LONG  WAY  TYPE  A  ZERO 
OTHERWISE  TYPE  AN  INTEGER 

0 

DO  YOU  WISH  TO  IMPOSE  MAGNITUDE  LIMITS? 

THE  DEFAULTS  ARE  0  =  MBRIGHT  .LE.  MASTEROID  .LE.  MFAINT  *  18.5 

TYPE  A  ONE  TO  ONLY  CHANGE  MBRIGHT 

TYPE  A  TWO  TO  ONLY  CHANGE  MFAINT 

TYPE  A  THREE  TO  CHANGE  BOTH  MBRIGHT  AND  MFAINT 

TYPE  A  FOUR  TO  CHANGE  NEITHER  MBRIGHT  NOR  MFAINT 

4 
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YOU  ENTERED  4 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

etc. 
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TABLE  A5 

ENTER  THE  MONTH  AND  DAY  OF  INTEREST  (.PROGRAM  ASTpT  3/80) 
AS  IN  8  17  FOR  AUGUST  SEVENTEENTH 

IT  IS  ASSUMED  THAT  THE  YEAR  IS  THE  CURRENT  ONE 


2 _ 5 

YOU  ENTERED  2  5 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

2 


OUTPUT  IS  FOR  2  5  80  AT  6PM,  12AM,  AND  6AM  LOCAL 
IF  YOU  WANT  QUICK  INPUT  ENTER  A  ZERO 

0 

ENTER  MB.MF , IL , IU , EL , EU .TAPE  UNIT  » 

15.5  18  0  14  0.23  1 

THE  MAGNITUDE  LIMITS  ARE  15.50  .LE.  MASTEROID  .LE.  18.00 
THE  INCLINATION  LIMITS  ARE  0.00  .LE.  INCL  .LE.  14.00 
THE  ECCENTRICITY  LIMITS  ARE  0.00  .LE.  ECC  .LE.  0.23 
TAPE  SHOULD  BE  MOUNTED  ON  MTl 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 

0 

ENTER  MB.MF, IL, IU,EL,EU,TAPE  UNIT  » 

15.5  18.  0.  14.  0.  0.23  2 
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THE  MAGNITUDE  LIMITS  ARE  15.50  .LE.  M ASTEROID  .LE.  18.00 
THE  INCLINATION  LIMITS  ARE  0.00  .LE.  INCL  .LE.  14.00 


THE  ECCENTRICITY  LIMITS  ARE  0.00  .LE.  ECC  .LE.  0.23 
TAPE  SHOULD  BE  MOUNTED  ON  MT2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
kk 

TRY  AGAIN  -  BAD  READ 
ENTER  MB.MF.IL. IU. EL, EU. TAPE  UNIT  « 

15.5  18.  0.  14.  0.  0.23  2 

THE  MAGNITUDE  LIMITS  ARE  15.50  .LE.  MASTEROID  .LE.  18.00 
THE  INCLINATION  LIMITS  ARE  0.00  .LE.  INCL  .LE.  14.00 
THE  ECCENTRICITY  LIMITS  ARE  0.00  .LE.  ECC  .LE.  0.23 
TAPE  SHOULD  BE  MOUNTED  ON  MT2 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 
3 

THE  OUTPUT  IS  AUTOMATICALLY  GIVEN  BY  INCREASING  ASTEROID  NUMBER 

YOU  MAY  HAVE  IT  SORTED  BY  INCREASING  RA,  MAG,  OR  ANGULAR  SPEED 

THESE  SORTS  REQUIRE  SOME  EXECUTION  TIME  AND  YOU  HAVE  EIGHT  CHOICES 

TYPE  A  ZERO  FOR  NO  ADDITIONAL  SORTING 

TYPE  A  ONE  FOR  AN  RA  SORT 

TYPE  A  TWO  FOR  A  MAG  SORT 

TYPE  A  THREE  FOR  AN  ANGULAR  SPEED  SORT 

TYPE  A  FOUR  FOR  AN  RA  AND  A  MAG  SORT 

TYPE  A  FIVE  FOR  AN  RA  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SIX  FOR  A  MAG  AND  AN  ANGULAR  SPEED  SORT 

TYPE  A  SEVEN  FOR  ALL  THREE  SORTS 


71 


6 


YOU  ENTERED  6 

IF  YOU  WISH  TO  CHANGE  THIS  ENTRY  TYPE  A  ZERO.  IF  NOT  THEN 
ANY  INTEGER 


YOU  DESIRE  AN  APPARENT  MAGNITUDE  AND  AN  ANGULAR  VELOCITY  SORT 
EXECUTION  BEGINS  NOW 
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system.  The  user's  Interaction  with  the  software,  complete  documentation  of  the  calculations,  a  user's 
guide,  sample  outputs,  and  the  rationale  of  observing  asteroids  are  also  discussed. 
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